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Abstract 

Quantum dynamica of a massless Dirac particle in time-dependent ID box and 
circular billiard with time-dependent radius is studied. An exact analytical wave 
functions and eigenvalues are obtained for the case of linear time-dependence of the 
boundary position. 

Neutrino billiard was discussed first by Berry and Mondragon in the 
context of time-reversal symmetry breaking and quantum chaos in relativis- 
tic systemsfl]. Recently 2D, massless Dirac equation has attracted much 
attention in the context of graphene physics. Graphene is a planar mono- 
layer of carbon atoms tightly packed into a two-dimensional honeycomb 
lattice. The quasiparticle motion in graphene in the low-energy regime is 
described by zero-mass Dirac equation with different (bulk and nanoribbon 
graphenes) boundary conditions [2, 4, 3]. 

Also, so-called graphene billiard has become attractive topic [5] in 
the context of graphene physics. 

In this work we treat the problem of neutrino billiard with mov- 
ing boundaries. Besides its relevance to graphene billiard such system 
is of importace in the context of relativistic quantum Fermi acceleration, 
too. Studying of time-dependent quantum requires solving of the two- 
dimensional Schrodinger or Dirac equations with time-dependent boundary 
conditions. We note that quantum mechanical wave equations with time- 
dependent boundary conditions have been subject of extensive study dur- 
ing past three decades [6] -[12]. The early treatment of the one-dimensional 
Schrodinger equation with non-stationary boundary conditions dates back 
to Doescher [6] who explored quantum dynamics of a particle in square 
well with moving walls. Later Munier[7] and Pinder [8] studied this prob- 
lem for some special cases. More comprehensive and mathematically strict 
treatment of the Schrodinger equation with time-dependent boundary con- 
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ditions has been done in a series of papers by Makowski et al [9]-[?]. 
Seba has studied this problem in the context of quantum Fermi accel- 
eration [11]. The case of time-periodic boundary condition is treated by 
Scheiniger and Kleber [12]. Three-dimensional case was treated by Yuce 
for the case of spherical boundaries [?] Despite the progress made in the 
study of Schrodinger equation with time-dependent boundary conditions, 
conditions, relativistic extension of this problem for the case of Dirac equa- 
tion with moving boundary conditions has not yet been considered. In this 
work we address the problem of massles Dirac particle in time-dependent 
confined geometries considering most simplest cases, time-dependent ID 
box and time-dependent circular billiard. For the case of special time- 
dependence of the boudnary conditions we obtain analytically eigenfucn- 
tions and eigenvalues of the Dirac equation for massless particle. 



1 Time-dependent one dimensional box 

Consider the following time-dependent Dirac system {% = c = 1) 

i-Zj: = (ap + rnfi)^, (1) 

with fy(t, x) = (^i(t, x), ^2(t, x)) T being two-component spinors, in the 
domain 

D = {(t, x), < t < T, < x < L(t)}, 
where the right boundary L(t), is time-dependent and 
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P = 

ox 

and 

are the Dirac matrices, for which the boundary conditions are given as 

V 1 {t,Q) = Q, V 1 {t,L(t)) = Q, 0<t<T (2) 

The Dirac equation for (fixed) boundary conditions, given in the one- 
dimensional box, ( L{t) = const) was treated in refs. [13, 14] where the 
solution of the Dirac equation 

{ap + mc 2 f3)^ = E^l 



2 



is obtained as 



~ E^+mc? C0s(/C n x) 

with being the normalization constant and k n = n7r/L,n = 1,2,..., 

E n =K + (mc 2 ) 2 ] 1/2 . 

To solve the eq. 1 we will restrict ourselves by considering the mass- 
less case, i.e., we assume that m = 0. The eq.l with the boundary con- 
ditions 2 cannot be solved until time-dependent boundary conditions are 
not replaced with fixed ones. To do this we use in the eq.l the following 
substitution: y = that reduced the domain D to D\ = {(t, x), < t < 
T,0<y< 1}. 

In this case the eq. (1) can be rewritten as 

) dt dy » dy ' o 

z^ = -L- 1 ^ + iLL- 1 ^ ' [ } 

where L = and the conudary conditions are given as 

*i(*,0) = 0, ^i(t,l)) = 0, (4) 

Time and coordinate variables in the eq.5 cannot be separated for 
arbitrary time-dependence of L(t). The only case for which variables can 
be separated is L{t) = at + b. 

Introducing new time variable 

ds 1 j fat + 6\ 



L(s) a 

and using the substitutions ^i(r, y) = e~ lXT f(y), ^i(r, = e~ tXT g(y) we 
get from the eq. (3) 



This system can be reduced to the second-order equation as 

d 2 _ , . , d 



(5) 



(i - «V)-^ - 2a y( lX + - iaX )f = -^f- ( 6 ) 

Using the substitution 



I -ay 
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the eq. (6) for \a\ < 1 can be reduced to the following equation 

z{l - z)f" + (c - (a + b + l)z)f - abf = 0. 

whose solutions are hypergeometric functions. 

In our case the hypergeomtric series can be written as 

F x {z) = F( ai , (3 h 7i, z) = £ { -^rZ k = (1 " zf 1 - 

Second solution can be found as 

F 2 (z) = A 

Then the general solution of the eq. (6) can be written as 

where A 1 = exp(iln(a)) 

Then from the boundary conditions given by the eq. (4) we get 

A - B ' (V) -(T) -°- < 8 > 

Solving this equation we obtain the eigenvalues as 

A n = 2airn | In |1 — a\ — In |1 + a\\~ l . (9) 
Corresponding eigenfunctions are 



U(y) = M ((l - m/p^ - (l + m/P^ 

= ^ ((1 - m/r A " /a+1 + (1 + ay)" iA "/ a+1 ) - (10) 
layM ((1 - m/p A " /a - (1 + ayy lK/a ) 

We note that for a — > the eigenvalues coincide with A n — > 7m, and 
eigenfucntions become 

fn(y) -> Asin{imy), 
9n(y) -> -^cos(irny). 

Then the complete set of solutions of the eq. (1) for linearly moving 
boundaries can be written as 



\f! ln (t,x) = Mexp 
V2n(t,x) = Mexp 



iXn l n ( at+b \ 
a \ b J 

iXn | n ( at±b\ 
a V b ) 



i^n/a / 1 ax \-iA„/a N 



/ 1 ax \ lA n/a / ax \ 

V 1 V 1 at+fey 1 



f 1 //-i _o^\- iA «/«+ 1 , /-i , _a^_\- iX n/a+l\ _ 
' \ iX n -a {K 1 at+b) T V 1 ai+&J / 

(ii) 

Thus we have treated the Dirac equation for massless particle in 
one-dimensional infinite square well with a time-dependent wall. 

2 Time-dependent circular billiard. 

Consider the Dirac equation with the boundary conditins given at the circle 
with time-dependent radius x 2 +y 2 < r^t). Radial equation can be written 

as 



(12) 



j SQ(T,t) _ dP(r.t) _ hp, 

dt dr r V ' /' 

where k is an integer number. Boundary condition is given by 

P(r (*),*) = 0. (13) 

Using substitution 



y = — r = / 



ds 



the eqs.(12) and (13)can be rewritten as 

i^ = «™^ + ^-*«(v,r), 



(14) 



where the boundary condition is time-independent now and given as 

P(M) = 0. (15) 

Separating time and coordinate variables in the eq.(14) for fo(t) = 
and ro (t) = at+b and reducing the obtaibned first-order system into second 
order equation we have 



2/ 2 («V " + 2a(a + i\)y 3 f- + [(ia\ - X 2 )y 2 + k(k + 1)]/ = (16) 

ay 1 ay 
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Solution of this equation can be written in terms of (regular at y = 
O)hypergeometric function as 



fx(y) = V 1 F (a, a + 7, a 2 y 2 ^j , 



where 



a 



k 



iX 1 

H h - 7 = 

2m 4 f 



+ 1. 



For secnd component we have 



g x (y) = (A - z(l - k)a)y k j - zay/ A (y) 

o s 



The eigenvalues X n can be found from the boundary condition: 

.2\ _ n (17) 



'♦i 



iX 1 1 



'4 



z'A 3 
+ 2^ + 4^° =°- 



2m 4' 2 

For k = the solution of the eq.17 can be found analytically: 

A n = 27m a/ in 

\1 — a/ 

The solutions of the eqs.(12) and (13) for k = can be written as 
P n (r,f) = N exp I - 



iX n (at + b\\ 



Q n {r,t) = Nexp 



iX r 



In 



'at + 6 s 



2a ""V 6 77 \ro(t) 
where N is the normalization constant. Thus we have obtained analyti- 
cally the solution of the massless Dirac equation for time-dependent ID- 
box and circular billiard. The circle is considered as monotonically ex- 
panding (contracting) with constant velocity. For general case (e.g., for 
breathing circle etc.) the solution can be obtained numerically. The above 
system is of importance because of its relevance to quantum Fermi accel- 
eration in relativistic systems and particle transport graphene as in real 
situation the boundaries on such system is not strictly fixed. Another 
important aspect of the above treated problem is its relevance to Dynam- 
ical (fermionic) Casimir effect. In this context the extension of the above 
problem to the case of non-integrable billiard geometries and time-periodic 
boundaries could be important. Currently such studies are in progress. 
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